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Abstract 



We reconsider the evolution equations for multigluon correlators derived in |l2j ]. We show how 
to derive these equations directly in terms of vector potentials (or colour field strength) avoiding 
the introduction of the concept of colour charge density in the intermediate steps. The two step 
procedure of deriving the evolution of the charge density correlators followed by the solution of 
classical equations for the vector potentials is shown to be consistent with direct derivation of 



evolution for vector potentials. In the process we correct some computational errors of 12 and 
present the corrected evolution equations which have a somewhat simpler appearance. 



1 



1 Introduction 



With the advent of HERA there has been a great increase in the scope of the theoretical effort 
to understand the physics of hadronic scattering at high energy. This is a challenging subject 
especially as it might provide a bridge between perturbative partonic physics of short distance 
processes (e.g. DIS at moderate XBj) and soft physics of hadronic states which, presumably, 
dominate the high energy asymptotics. The borderline between the two — the "semihard" physics 
— is interesting also in its own right. In essence it is the physics of partonic systems which are, 
on one hand dense enough for new collective phenomena to play important role but, on the other 
hand are perturbative since the average momentum transfer between the partons is high enough. 
In this semihard regime one expects to see perturbativelly controllable nonlinear effects which 
depart from the standard linear evolution of DGLAP or BFKL Q type and subsequently lead 
to unitarization of hadronic cross sections. 

An approach to this high partonic density regime from the partonic side has been spearheaded by 
Levin and collaborators H based on earlier work by Mueller [|J and is an " all twist" generalization 
of the GLR recombination picture ^,0] . It lead to the formulation of a nonlinear evolution equation 
which exhibits a perturbative mechanism of unitarization. Analysis of this equation suggests that 
already at present HERA energies the nonlinearities in the gluon sector are considerable and linear 
evolution for gluons should break down. Better experimental data on gluon distributions would 
be extremely valuable in order to verify/falsify this assessment. 

A complementary approach was pioneered some years ago by McLerran and Venugopalan It 
was later somewhat reshaped conceptually and considerably developed technically in a series of 
papers [|l0- 1 3 1 . Here the idea is that in the high density regime rather than using partonic language, 



it is more appropriate to use the language of classical fields. The hadron then is considered as an 
ensemble of configurations of the gluon field. The statistical weight that governs the contributions 
of different configurations to the ensemble averaging changes when one probes the hadron on 
different time scales. Decreasing xsj corresponds to increasing the time resolution and therefore 
corresponds to probing the hadron on shorter time scales. The change of the statistical weight 
with x Bj is governed by the evolution equation. As long as the field intensity is large enough this 
evolution should be perturbative in a s but essentially nonlinear in the field intensity itself. 

This evolution equation to first order in a s was derived in |ll] , [l2| . We will refer to it as the JKLW 
equation in the following. In []l3|] the double logarithmic limit of this evolution was considered. It 
was shown that in this limit the evolution of the gluon density becomes unitary at large density. 
Qualitatively, the evolution is very similar to that discussed in || , although the details are different. 
A detailed numerical study of the doubly logarithmic limit of the JKLW evolution was recently 
performed in fl4|| . 

Technically the derivation of [^2[ is fairly involved. Several consistency checks were performed in 
jnj and |l3| to make sure that the known results are recovered from the general evolution equation 
in the weak field limit. This includes the BFKL equation, the doubly logarithmic approximation 
to the DGLAP equation and the GLR equation. It is, however, desirable to have some additional 
independent checks on the equation which do not involve the weak field limit. It is the aim of the 
present paper to provide one such check. 

In a nutshell the issue we address is the following. The evolution equation in |Tl] [l| was derived 
invoking a two step procedure. Rather than considering directly the evolution of the correlators 
of the gluon field, one first considers the evolution of the colour charge density. In the second 
step, one re-expresses the evolution equation for the charge density as the evolution equation for 
the vector potential (chromoelectric field). The evolution of the field correlators is in fact what 
one is after, since it is the vector potential and not the charge density that couples directly to 
fermions and, therefore, are more directly related to physical observables. Our observation in this 
paper is that one can avoid the introduction of the colour charge density altogether and derive 
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the evolution equations directly for the field correlators. This procedure has the advantage of 
being somewhat simpler both technically and conceptually. Nevertheless, the final result for the 
evolution should be the same as in the two step procedure of 10 The comparison of our 



results with the earlier derived formulae provides us with a consistency check on the calculation. 

We find, in fact, that the results presented in |l2| are not entirely correct. However, after cor- 
recting some algebraic mistakes in Jl2| we show explicitly that the two approaches yield identical 
results. We provide the corrected expressions for the "kernels" of the evolution equation, which are 
somewhat simpler than the expressions found in jl2| . We also clarify the issue of possible Gribov 
ambiguity and show explicitly that the divergent Jacobians, which appeared in the intermediate 
steps of the derivations in p3| , cancel completely in the final expressions for the correlators of 
the chromoelectric field. Therefore, the Gribov ambiguity, although affecting the relation between 
the colour charge density and the chromoectric field, does not affect the evolution of the field 
correlators, at least to order a s . Since the procedure discussed in the present paper avoids the 
introduction of the charge density entirely, the whole approach is free from the Gribov problem. 

Perhaps somewhat surprisingly the corrections to the results of that we find do not affect 
either the weak field limit discussed in jllj or the doubly logarithmic limit of |lj| . They therefore 
have no bearing on the derivation of the BFKL equation in our approach and also do not help 
to reconcile the doubly logarithmic limit of the JKLW equation |lj| with the nonlinear equation 
studied in aft. 

The plan of the paper is the following. In Section 2 we briefly recap the procedure of the derivation 
of the evolution equation as described in [Il|-|l3(] and reformulate it directly in terms of the gluon 
field correlators. In Section 3, using some of the results of [Q, we calculate the real and virtual 
parts of the evolution in terms of the field correlators and provide the corrections to the results 
of p^ |. Finally, in Section 4 we discuss our results. 



2 The JKLW evolution 

First, let us briefly recall the framework and the results of In this approach the averages 

of gluonic observables in a hadron are calculated via the following path integral 

< 0(A) >= J DpDA fi O(A)exp{- J d 2 x x F[p a {x x )] - i J dS-trF"" F M „ 

+ jj- J d 2 x ± dx~5(x~)p a (x±)trT a W- 00:00 {A~}(x~ ,x±)} (1) 

where the gluon field strenght tensor is given by 

Fr = d»A v a - d»A» - gf ahc A»A v c (2) 
and W is the Wilson line in the adjoint representation along the x + axis 



W-oo,+oo[A }(x ,x±) = Pexp 



ig J dx + A a (x + ,x ,x±_)T a 



(3) 



The hadron is represented by an ensemble of colour charges localized in the plane x~ — with 
the (integrated across x~) colour charge density p(x±). The statistical weight of a configuration 
p(x±) is 

Z = cxp{-F[p}} (4) 
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In the tree level approximation (in the light cone gauge A + = 0) the chromoelectric field is 
determined by the colour charge density through the equations 

F +l = -8{x-)a\ Xl _) (5) 
9 

and the two dimensional vector potential a l (x±) is "pure gauge", related to the colour charge 
density by 

d l ai - d j a\ - f abc a\a? c = 

d l al = - Pa (6) 



Integrating out the high longitudinal momentum modes of the vector potential generates the 
renormalization group equation, which has the form of the evolution equation for the statistical 
weight Z0Q 

— Z = a s \ — - , r r . . \Zy{u, v)] — - 7 . (7) 
d( \ 2 5p(u)Sp(v) 1 y n Sp(u) 1 n ] K ' 

In the compact notation used in Eq. (^) , both u and v stand for pairs of colour index and transverse 
coordinates, with summation and integration over repeated occurrences implied. The evolution in 
this equation is with respect to the rapidity £, related to the Feynman x by 

C = lnl/£B (8) 

Technically it arises as a variation of Z with the cutoff imposed on the longitudinal momentum 
of the fields A 11 . The quantities \[p] an d cr[p] have the meaning of the mean fluctuation and the 
average value of the extra charge density induced by the high longitudinal momentum modes of 
. They are functionals of the external charge density p. The explicit expressions have been 
given m 12 and it is our aim in this paper to provide a check on these expressions. 

Eq. (Q) can be written directly as an evolution equation for the correlators of the charge density. 
Multiplying Eq. (Q) by p{x\)...p{x n ) and integrating over p yields 

< p(xi)...p{x n ) > = 

< p(xi)...p(x m - 1 )p(x m+ i)...p(xk-i)p{xk+i)...p(x n )x(x m ,Xk) > 



0<m<fc<n+l 



^ < p{xi)...p{xi-x)p{xi + x)...p{x n )a{xi) > 



0</<n+l 



(9) 



This set of equations for the correlators of the colour charge density completely specifies the 
evolution of the hadronic ensemble as one moves to higher energies (or lower values of x). 

The evolution equations for the correlators of the charge density can be rewritten as equations for 

1 A11 the functions in the rest of this paper depend only on the transverse coordinates. For simplicity of notation 
we drop the subscript _L in the following. 
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the correlators of the vector potential |l3| . 



^ < <(zi)-<>n) > = 



< <(^0---4tU^-i)<^(^+i)"-<(^)<(^) > 

0<Kn+l 

+ <<^)"<: 1 1 (^-i)C + 1 1 (wi)- 



0<m<k<n+l 



X <\7_\(xk-i)a^ + \(x k+1 )...all(x n )xlz%( x ™,Xk) > 



(10) 



The quantities \ab an d a a have a very simple physical meaning. The high momentum modes of 
the vector field which have been integrated out in order to arrive at the evolution equation induce 
extra colour charge density dp. The average value of this induced density and its mean fluctuation 
appear in the evolution equations eq.(^|) as er Q and x a b- The appearance of the induced colour 
charge density leads to the change in the value of the chromoelectric field through the solution of 
Eq. @ with p + dp on the right hand side. The quantities a l a and \^ b are the average value and 
the mean fluctuation of the induced field respectively. 

It is perhaps helpful to explain how a\ and xL> were obtained in As shown in |p^| , the 

induced charge density can be decomposed into two pieces ^] 

5p = 6pi + 5p 2 (11) 

The first piece 5pi is order g while the second piece 8p 2 is order g 2 . The 5p\ is time dependent, and 
has zero average value, while its mean fluctuation is order g 2 . The Sf>2 being 0(g 2 ) contributes only 
to the average value of Sp and not to the mean fluctuation. Assuming that the classical equations 
eq. (|^) hold not only for the background field but also for the relevant part of the fluctuation field 
one can solve those equations perturbatively. Writing 

6a z = 5a\ + 6a 2 (12) 

with Sai being 0(g) and Sa 2 being 0(g 2 ) and keeping in the classical equations all terms to order 
g 2 we have 

D l ab Sa{ b - D J ab 5a\ b + D\ b 5a{ b - D J ab 5a 2b - f abc Sa\ b 5a\ c = 

d l Sa\ a + d l 8a\ a = -(Sp la + 6p 2a ) (13) 



We have defined for convenience 



To order g we find 



"afc = fabc^l 

Dl b = d l 6 ab + al b (14) 



Sa\ = -D i ±8p 1 (15) 



Therefore, to order g 2 

Xab( x ' V) = r lc(x, u)Xcd(u, v)rlj> b (v, y) (16) 



The reason for the notation p rather than simply p will be explained in the next section. 
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with 

ri b (x,y) = -<x\[D i -^ B ] ab \y> (17) 

Here < x|C% > denotes a configuration space matrix element in the usual sense. 
At order g 2 we have 

5o ^ = ~ Dl j5 6h ~ l eijdi id Sai x Sai (18) 

Here the cross product is defined as 

A x B = fabc^A\B{ 

We thus have 

^afc) = r ab(^ + PobcC 35 ! U > w )X6c( u , «) (19) 

with 

pi bc (x, V> z) = -^dn^} ad )(x,u)f dfe e kl r k fb (u >y y ec (u,z) (20) 

The procedure of deriving eq.(|i"o|) employed in jl2],[l3| consists, therefore, of two steps. One first 
splits the gluon field into the classical background field a' 1 and the fluctuation field a' 1 . The modes 
of the fluctuation field with longitudinal momenta in some range a s In -^-r are assumed to be small. 
One defines operatorialy the induced charge density Sp in terms of the fluctuation fields and 
the quantities a and x are calculated by integrating out the fluctuation fields perturbatively. In 
the second step, one solves classical equations of motion which include the induced charge density 
and calculates a % and \ % K 

Clearly, consistency requires that the two step procedure that leads from eq.(|l]) through eqs.(|7|||) 
to the evolution equations eq.(|l(], [16], |l^) must be equivalent to the following. Start with the 
equivalent of eq. (|l|) 

< 0(A) >= J Da t DA' i O(A)Z[a l (x ± )}exp{-i J d^x^trF^F^ 

-■fir J d 2 x ± dx'6(x-)d t al{x ± )tTT a W. 00 ^[A-](x- 1 x ± )} (21) 

Integrate out the high longitudinal momentum components of a M as before, but instead of calcu- 
lating the induced charge density a a and \ a b^ calculate directly the induced chromoclectric field 
a l a and Xab- Technically this calculation is somewhat simpler, since there is no need to consider 
the operator Sp, which is nonlinear in the fluctuation field a M . Instead, one directly calculates the 
distribution of the static component of a M . The resulting evolution equations should coincide with 
eq-©- 

With this formulation one circumvents completely the need to introduce the colour charge density 
p and to solve classical equations for a 1 in terms of p. While one may want to introduce p for 
reasons of convenience, our present understanding is that it is not necessary from the point of 
view of physics. The physics that our approach is meant to address is that of the evolution of the 
hadronic ensemble. The relation between a 1 and p on the other hand is supposed to hold at every 

3 We note that the statistical weight Z[ai\ which appears in eq. (|2l|) is not equal to Z of cq.(^|) since going 
from eq. (hi) to eq. (fell) involves the change of variables p — » a;. The two statistical weights, therefore, differ by an 
appropriate Jacobian. 
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value of £, and therefore itself is unrelated to evolution in £. The concept of p may be sometimes 
useful to formulate models for the statistical weight Z at some particular value of Q as was the 
original motivation of Q. This could then serve as an initial condition for the evolution. This, 
however, is a separate question and we do not intend to address it here. 

Before we proceed further, we wish to make one more comment about the relation between the 
chromoelectric field and the colour charge density eqs.(|l5|,[l8|). Both these equations contain the 
dangerous factor (dD)^ 1 . The operator dD has zero, as well as negative eigenvalues and is very 
reminiscent of the operators usually associated with the Gribov ambiguity in nonabelian gauge 
theories. In fact, it is quite clear that it has precisely the same origin. The second equation in 
eq.(|l^) has the form of the Lorentz like gauge fixing condition on the fluctuation field 8a. Since 
the calculation is performed in a nonvanishing background field, the Lorentz gauge indeed suffers 
from Gribov ambiguity precisely due to negative eigenvalues of the operator dD. 

Given this, one may worry that our pcrturbative calculation is plagued with the Gribov ambiguity 
^. However this is not necessarily the case. The point is that 8p itself is not arbitrary. It is 
calculated through the fluctuation field and, at the end of the day, is averaged over with some 
statistical weight Z[Sp\. It could well be that the statistical weight is such that it only allows 
induced charge density of the form Sp = dDX with regular X . If that is the case, the dangerous 
denominator cancels and the induced field is well defined and regular. In fact, in our present 
formulation where the calculation is performed directly in terms of the field, it is almost clear 
that this should indeed happen. In this setup one calculates directly 8a, and eqs.(15 Wf) should 



be read from right to left, as equations determining an auxilliary quantity 5p through 5a rather 
than the other way round. The operator dD then appears in the numerator and all expressions 
are regular. In fact we will show in the next section by explicit calculation that all "dangerous 
factors" indeed cancel in the final expressions for \lb anc ^ °a- 

Note that, if one insists on formulating the problem in terms of the colour charge density, the 
absence of the Gribov ambiguity implies a nontrivial consistency condition on the statistical weight 
Z[p\. Taking an arbitrary weight Z will render the calculation of chromoelectric field correlators 
ill defined especially at strong fields (strong coupling). This was indeed observed in the numerical 
calculation fl5| where a simple Gaussian in p was used as the weight function 0. 

In the next section we will calculate a l a and \ab induced by high longitudinal momentum modes. 



3 The Induced Chromoelectric Field 

The main ingredients needed for the calculation of the induced chromoelectric field are the eigen- 
functions of the quadratic action for the small fluctuations in the static background a 1 . Solving 



the classical equations of motion that follow from the action eq.(21) at fixed a 1 we find (in the 
gauge d l A l {x + — » — oo) = 0) the classical solution 



A d =0, 4, = a^MO (22) 
Defining the quantum fluctuation field a M by — A^ + a M and expanding the action to second 

4 In standard perturbation theory, the Gribov ambiguity does not show up in any finite order. This is due to 
the fact that one expands the operator dD and its inverse in powers of the coupling constant. To leading order 
then the operator does not have any negative eigenvalues, which ensures that no problems arise in finite order 
perturbative calculations. Our situation is, however, different. Since our background field is not assumed to be 
O(g), the operator cannot be expanded. Therefore, there is no guarantee that the problem does not show up even 
in perturbation theory. 

5 This problem does not arise in the more recent numerical work [h™ since in effect this work uses a different 
definition of p for which the relations analogous to eqs.(|15l[ia) do not involve singular factors. 
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order in a M we have 



S = —jl a-K xy a- + 2a-(d+Da - 2 fa) + 2d + a i d~ a 1 + a 1 



D-d" - IV l)> 



Here we are using the notation 

[fa] a (x + ,x~,Xj_) = S(x~)al b (xx)al(x + ,x~,x±) 
Da = L>>K = {d l S ab + e{x-)o^ ah )a\ 



and as previously 



r I 

Jab c a c 



(23) 



(24) 
(25) 



The operator K is 



K xy 

^ab 



{d+Y5 ab + d l a\ b 8{x-) 



d- 



(26) 



Note that there is no ambiguity in the definition of the operator l/d~ in this expression. It is 
defined in the sense of principal value. This follows directly from the fact that the matrix a l ab is 
antisymmetric and therefore the term involving l/d~ in eq.(^3|) vanishes for zero frequency fields. 

This eigenfunctions of the quadratic fluctuation operator have been found in fl^| and we cite here 
the relevant results. The calculation is performed in the lightcone gauge A + = with the residual 
gauge freedom fixed by the condition 



It is convenient to define an auxiliary field 



-co) = 



or = a- + K- l {d + Da - 2 fa) 
This field can be seen to decouple from <Zj. Its correlator is 

< K a~ >= K~l 

The operator K Eq.(p6|) has zero modes. Defining the projector matrices 77 and \i by 



and 



/i + 7] = 1, /i = fX, T] = 1] 

we can write the normalizable zero modes of K in the form 

fa(x±,X~ ,p~) = {l ab f(x ± ,p~) 



(27) 

(28) 
(29) 

(30) 

(31) 
(32) 



The operator K is therefore, strictly speaking, non invertible. The operator K^ 1 in eq.(|28|) has to 
be understood as the inverse of K on the space of functions which does not include the functions 
Eq. (|32"|) . Further, it is only the nonzero mode part of a - that enters the definition of a~ in eq.(28). 



For our calculation we will need the properly normalized solutions of the equations of motion that 
follow from the action eq. (E3|) . The complete set of these solutions was found in |12f 



p ? r 



ge" 



d 2 pi 



i cxp ( i^Z x ~ ip±x± ) r[ , (p ) 



+ 9(x~)U(x ± ) 



) cxp f - ip±x±) [U^v+A (Px) + 9(x )Y+, r 



(33) 
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The frequency p~ is a good quantum number since the background field is static. Here r is the 
degeneracy label, which labels independent solutions with the frequency jT . In the free case it is 
conventionally chosen as the transverse momentum, {r} — {p 1 }. The matrix U(x±) is the SU(N) 
matrix that parametrizes the two dimensional "pure gauge" vector potential a l (x±) 

a'CsO = -U(x x )d i U^(x x ) 
9 

The auxiliary functions 7^ , v± are all determined in terms of one vector function. Choosing this 
independent function as v l _ we have 



J +,r 



r+,- 



2D 



(34) 
(35) 



where we have defined the projection operators 

/;'/)' 

t ij = 5 i3 



D 2 
gigj 



d 2 



L lJ = 
ft = 



D l D 3 

QiQ3 



(36) 



The proper normalization of the eigenfunctions requires v l _ to be chosen as complete set of eigen- 
functions of the two dimensional Hermitian operator O^ 1 



[(t-l)0- 1 (t-l)Z(x ± ,y ± ) = 

<^l€,-2 



d 2 



D 21 



D 2 



\y± > 



(37) 



such that 



d 2 r±v l _^ a (x ± )v*_ J rb (y ± ) = 



[0-^ b (x x ,yx) 



(38) 



The rotational scalar operator S is 



1 d i D i d l D l 
q _ _ 1 of- — - — 1 

D 2 [ d 2 D 2l[ d 2 D 21 



(39) 



For further use we also need the expression for the a component of the fluctuation field. Using 
the explicit expression for the operator K from p2[ we get from eq.(Eq) 



a (x ,x±,p ) = a —6{x ) / dy D l (a l -^ l + ) 



-e(-x-) 



dy-D\a l / dy~d l a 



2ip-v[^-^](t-l) ij vL(x ± ) (40) 



We note that this expression differs by a a; -independent constant from the one given in |l2| . The 
reason is that in fl2|| a constant have been subtracted from a~ such that f^°° dx~a~(x~ ) = 0. 
This corresponds to the symmetric definition of the integral in eq. (|io|) . This is incorrect, since it 
violates the residual gauge fixing d l a l (x~ — > —00) at the one loop level. We will see this explicitly 
later in this section. At any rate, straightforward albeit somewhat tedious calculation gives the 
result eq.(^0|) and this is the expression that will be used in the rest of this paper. 
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So far the formulae presented in this section (except for the corrected expression for a~ eg. (f40|) ) 
are identical to those that appear in jl2) with the only difference that the background charge 
density p has been substituted by the background field via p — —d % a % . Now however we will take 
a different route. Our aim is to calculate the order 0(a s lnl/a;) correction to the background 
chromoelectric field eq.(|^) directly, rather than to the background charge density. According to 
the discussion in the previous section (see also |l^]), we are therefore interested in the following 
two quantities 

a s lnl/xXab( x -l,yx) =< a % a (x±,x~ -> oo, x + )a 3 b (y±, y~ -> oo,x + ) > (41) 
a s Inl/x <r l a (x± ) =< a l a (x±,x~ — > oo,x + ) > (42) 

ft should be noted that, since the background is static, none of the quantities defined above depend 
on x + . 



3.1 The real part - the mean fluctuation 

It is a straightforward matter to calculate \lb- R- ecau that we need this quantity to order g 2 . 
The fluctuation fields a M are formally of order g themselves, and therefore to calculate the mean 
fluctuation we do not have to include loop corrections. Examining the expression for the general 
solution eq.(|33|) we see that it contains oscillating pieces, which do not contribute to the value 
of the field at infinity as well as the 7+ piece, which does not vanish at infinity and, therefore, 
determines the distribution of the vector potential there. 



Xab(x±,V±) =4tt / dp <j l +Jx±,p )7+, 6 (y±,-p )> 



(43) 



Using the explicit expressions for 7I we find after some trivial algebra 



Xl(x±,y±) = 2 < x ± \{^{D 2 - S- 1 ] — } ab \y x > 



(44) 



We now want to compare this with the corresponding result of |12| . The induced charge density 
dp in [pf is 



dp = Spi + Sp 2 



(45) 



with 



5pia(xj 



1 



ai( x - = 0) - / dx~d + al(x~) 
Jo 



fated* al( Xl _) / dy+ 



6(y + - x + ) - 6(x + - y + ) 



a-(y+,x ± ,x- =0) (46) 



and 



Sp2a{xj_) = fabc I dx [d + a\ {x)]a\ (x) 

1 



(9 l a^(x^) dy + a c (y + ,x±,x =0) dz + a d (z + , x±, x = 0) 



facefbdeO(z+ - X+)6(x + - y+) + f ab efcde0{x + - Z+)9(z+ - y+) 



(47) 
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(a) (b) (c) 

Figure 1: One loop tadpole diagrams contributing to a l a . The tadpole is calculated at x~ — > oo. 
Only Spi contributes to \. Substituting the expressions for a % and a~ into eq.(Hq) we find 



Thus, we obviously have 



S Pl = -2(dD)[^--^}(t-l)v- (48) 



7+ = --D*^Pi ( 49 ) 



This reproduces exactly eq. (|l5|) . Obviously the relation between x a b an< ^ Xa6> ec f-(0) i s a l so 
reproduced by this result. 

We note that our result for x a b 1Sl somewhat different than the one presented in p2[ . As discussed 
before this is due to an incorrect treatment of the ^"-independent component of a~ in p2|| . 

3.2 The virtual part - the average value of the field 

We now proceed to calculate the virtual part of the evolution kernel. For this purpose we have to 
calculate the zero frequency part of the {ij} and {i— } components of the fluctuation propagator. 
The calculation of the {ij} at zero frequency is straightforward. The result is 

lim < 3 t ab( X ~iy~' X -L'V-L>P~) = lim < a a( x ~> x ±iP~) a b(y~>y±>p~) > = 

p -^0 p~ — >0 

- i6 ij S(x~ - y-) 9(x~) < x ± \(-^) ab \ y± > +d(-x~) < x ± \(-^) ab \y ± > (50) 

The {i— } component is then calculated immediately using this result, eq.(|40|) and the fact noted 
earlier that the field a decouples from a 1 . The result is 

>= 



lim G lb( x 'V ; x ±>y±,p )= lim < a a( x >%±,p )%{y ,vx,p ) 

p~^0 p- — »0 

{j^)ab\y± > +0{-X )<-_L|V^2 



r D i _ d i 1 

i0(x -y ) 9(x ) < x JL \(— I ) ab \y ± > +9(-x ) < x ± \(-^) ab \y± > (51) 



We are now ready to calculate a l a . It is given by the one loop tadpole diagrams of Fig. 1. The 
vertex lc comes from the expansion of the Wilson line term in the action to third order in the 
fluctuation. The separate contributions of the diagrams can be written in terms of the fluctuation 
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propagator G Miy =< a^a" > in the following form 



la= - / dy-d 2 y ± G%(x-,y-,x ± ,y ± ,p- = 0)ej k Dlf^ e e mn G^(y-,y- ]y± ,y ± ,y+,y+) 
lb = -i I dy-d 2 y x G^(x- ,y~,x x ,yx,P~ = Q)fbcddp =y -G } c J d {y~, y~; y ± , y±, y + , y + ) 



lc = jy- I dy d A y±dy + dw + dz + 



x 8{y )(d l al(y±_))G l ac (x + ,y + iX ,y ,x_L,y±)G de (w+,z + ,y ,y ,y±,y±) 
(z+ - y+)%+ - w + )f bef f cdf - %+ - z+)9(z+ - w+)f bcf f def 



(52) 



The diagram Fig. la corresponds directly to the second term in eq.([T^). For this diagram we 
immediately find 



D 2 



{x,y)f bcd e kl Xcd(y,y) 



(53) 



ah 



The diagramms Fig. lb and lc correspond to the first term in eq. ( |l9| ) and can be written as 



D> 



<K(2)0) = < $P2 > 



(54) 



with 6 pi (cf eq.(|4 



' <V>2 •„= /„/„ / (d + al(x))a l c (x) > + ^ {facefbde ~ "abe j "cde) ' Oi\{x j_ ) 

x / r— ^ . dp" . 1 ,„ < a~(p~,x±,x~ = 0)a~,(-p~,x±,x~ = 0) > (55) 
A + te (p y 



Using the results for the equal time propagators from Jl2t] we obtain 



< Sp 2 >a = (facefbde ~ ^fabefcde) d i a\(x±) 



D 2 D 2 



1 a 
Da d- 2+ l 



+ fabc < X± 



+ R a (x ± ) 



t ij - l ij - 2a l d j 



^ - 2 (^2 -^2) S_ (^2 -^2) 



1 



T fe 4 _ _ 2 —d k c 



\X± > cd 



\X± >bc 



(56) 



with 

R a (x±) = /abc / d 2 y ± d 2 z ± 



2„,. j2„. d 2 pj_d 2 kj_ p\ ip JL ( XJ __ VJ j_ife JL ( a ; JU _z JL ) 



(2tt)* p2_ fc 2 

{ < VJ.I — ^ 

-C/(^i) <y±|C/ t [T*-L«] 



t ki - l kt ] \z± > 



^ ' 2{ & 2 '~D^ )S ^ { & 2 "D 2 ' 



[T kl -L kl ] U\z±>Ul(xx)\ 

J be 

(57) 
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Here the singularity in the integrand at p 2 ^ = k 2 ^ has to be understood in the sense of the principal 
value 

1 p\ - k\ 

p\-k\ ~ (p\ - k\ ) 2 + e 2 ' 

Our final result for the induced field is given by the sum of eq. ( f53| ) and eq. ( p4| ) (supplemented by 
eqs.pg)). 



4 Conclusions 

To summarize, the final results of this paper are eqs.([44|) and ( ]54| , |56| ). They supercede the corre- 
sponding results of JlJ] and (lj| . 

We now want to comment on this result. The first thing to observe is that the dangerous denom- 
inator 3D does not appear in these expressions. The Gribov problem mentioned earlier therefore 
does not affect our calculation, at least to order a s . 

The result for the induced field differs from the corresponding formulae in |[2| and in two 
ways. One reason is the improved treatment of a~ relative to p2|| . Now we are in the position to 
understand why the expression for a~ used in is inconsistent with the residual gauge fixing. 
In the previous section we have calculated the induced vector potential far at infinity x~ — > oo. 
It is not much more difficult to calculate it everywhere in space. Diagrammatically it is given by 
the same diagrams as Fig.l except the coordinate on the free end of the propagator is some finite 
x~. The difference in the analytic expressions eq.(|54|) is that the surface charge density 5p2 is 
substituted by the local charge density integrated up to the longitudinal coordinte x~ 



- e ^)^i j d v~ < s ti(y~) > + - °(- x ~)q2 J 



dy- < Sj+(y-) > (58) 



This expression makes it explicit that the induced field vanishes at x~ —* — oo. Therefore, our 
calculation clearly preserves the residual gauge condition d t a t (x~ — » — oo) = 0. However if we 
were to subtract the zero momentum piece from the field a~ as done in Jl2[ , the integration 
limits in eq.([io|) would become symmetric — > \ + J x °°). The effect of this would be that 
G l ~ (x~ ,y~) would not vanish at x~ — > — oo. It is then obvious that we would have d l a l {x" — > 
— oo) ^ 0. The expression obtained in the present paper does not suffer from this problem. It 
is consistent with the perturbative ie prescription for regulating the 1 /p + gauge pole used in the 
earlier work [ fl7|| . 

Another difference between our present result and |l2| is the appearance of D 2 rather than 3D 
and Dd in the denominators in eqs. (|5^ , |5^ ). This deserves an explanation. This is also related to 
another point we want to address. Comparing eq.([54|) with eq.([lS|) one could wonder whether the 
present method of calculation of is consistent with the two step procedure of |lj, ^3). It may 
look like the relation between the induced field and the induced charge density we obtained here 
(eq.(|54|)) is different from the equation eq.^) which was used in the previous work. This however 
is not the case. The reason is that the 0(g 2 ) induced charge density 8p2 which appears in eq.([T9|) 
is not quite the same as < 5p 2 > in eq. (|54|) . The 5p2 was defined as complete 0(g 2 ) contribution 
to the average of induced density. In other words 



6p 2 =< Spi + 6p 2 > (59) 

with <5pi ; 2 defined in cqs.( pti| , f47| ). As we discussed above, the fluctuating part of the operator Spi 
is of O(g) and therefore indeed 8p\ can be identified with 8p\. However, the vacuum average of 
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8 pi is 0(g 2 ) and does contribute in eg. (|59|) . It can be shown that 

< 5p u >= f abc a\ < a} c {x + -> oo) > (60) 

This extra contribution turns d l 8ot l 2 into D l 5a 2 in the second equation in eq.(|l3|) if we use 5p 2 
rather than 8p2 in its right hand side. Taking account of this we see that the procedure described 
in section 2 is consistent with eqs.( p|j5^ , |5^ ). 

In fjjl it was assumed that < 8pi >= and thus the extra contribution of eq.(|60|) was overlooked. 
This lead to an apparent noncancellation of spurious factors which as we see now, do indeed 
cancel in the final result. 

Importantly, the corrections we find vanish in the limit of weak field considered in |p"l|| and also 
in the double logarithmic limit, where the field is considered not necessarily weak but slowly 
varying in the transverse plane |fl3|. This can be seen in the following way. Comparing eq.(|40|) 
to the appropriate expression in ||l2(| we find that the difference between the two is proportional 
to p. In the weak field limit one only needs to know a~ to order 1 and therefore the correction is 
unimportant. For slowly varying fields all terms proportional to p are also negligible. Therefore, 
the real part - x IJ - hi these two limits is insensitive to the correction we found here. The virtual 
part - a 1 - does not contribute at all in the DLA limit. In the weak field limit the correction is 
negligible since < Spi >= 0((a 4 ) 2 ) and one only needs 8p to order ct . 
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